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The phase structure of a single self-interacting homopolymer chain is investigated in terms of 
a universal theoretical model, designed to describe the chain in the infrared limit of slow spatial 
variations. The effects of chirality are studied and compared with the influence of a short-range 
attractive interaction between monomers, at various ambient temperature values. In the high tem¬ 
perature limit the homopolymer chain is in the self-avoiding random walk phase. At very low 
temperatures two different phases are possible: When short-range attractive interactions dominate 
over chirality, the chain collapses into a space-filling conformation. But when the attractive interac¬ 
tions weaken, there is a low temperature unfolding transition and the chain becomes like a straight 
rod. Between the high temperature and low temperature limits, several intermediate states are ob¬ 
served, including ^-regime and pseudogap state, which is a novel form of phase state in the context 
of polymer chains. Applications to polymers and proteins, in particular collagen, are suggested. 


I. INTRODUCTION 

A linear homopolymer is made of a single type of a 
repeat unit. An important example is polyacetylene, an 
organic conductive polymer which is the paradigm mate¬ 
rial for fractional fermion number mm- Additional ex¬ 
amples, among many others, are poly-L-lysine and poly- 
L-glutamic. The former is a food preservative with po¬ 
tential for wider, even pharmaceutically relevant antimi¬ 
crobial effects [3] while the latter is used for drug delivery 
against cancer [3]. 

From a theoretical point of view, the concept of a 
homopolymer chain is a useful coarse grained approx¬ 
imation, even in the case of a heteropolymer that ex¬ 
hibits only approximatively repeating patterns: For a 
sufficiently long chain the distinct monomers are simply 
combined into appropriate subunits, to dispose of the in¬ 
homogeneities in the monomer species. For example col¬ 
lagen, which is the most abundant protein in mammals, 
displays a repeated glycine-proline-X pattern, where X 
is any amino acid other than glycine and proline. DNA, 
RNA and the Ca backbone of a protein chain are addi¬ 
tional examples where a homopolymer approximation is 
occasionally profitably introduced mm- 

Here the phase structure of a chiral linear homopoly¬ 
mer is investigated, in terms of a universal energy func¬ 
tion [EHni; a chiral polymer is one where parity is broken, 
the mirror image of a stable chiral polymer conformation 
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is in general not stable. A chiral polymer often has a 
tendency to form helical structures. For example in the 
case of proteins right-handed helical structures are more 
common than left-handed ones. There are also chiral 
proteins that can form different right-handed and left- 
handed structures. An important example is collagen 
for which the right-handed polyproline I conformation is 
more compact than the left-handed polyproline II. 

It is found that the phase structure of a chiral ho¬ 
mopolymer is more complex than that of a non-chiral 
one. In particular, a chiral homopolymer can be in a 
pseudogap state. Here we give a brief description of this 
phase, following [I4HI8] for some abstract statistical sys¬ 
tem which displays a continuous symmetry that is asso¬ 
ciated with a complex order parameter 

( 1 ) 


Here p and p are just modulus and phase of this ab¬ 
stract order parameter. Such an order parameter is of¬ 
ten present e.g. in models of superconductivity. Further 
in chapter HE (see eq. (30l) we give description of this 
complex order parameter in terms of polymer degrees 
of freedom. The restoration of the continuous symme¬ 
try commonly takes place so that the free energy of the 
symmetry breaking state with non-vanishing condensate 
p 7^ 0 becomes larger than the free energy of the sym¬ 
metric state where p vanishes. But the symmetry can 
also be restored by phase decoherence, even when p 7^ 0. 
This occurs when the phase (p in the order parameter be¬ 
comes disordered so that (e*'^) = 0. This implies that 
the expectation value of the order parameter also van¬ 
ishes, (T') = 0. The system is then in the pseudogap 
state pAUl8] . The pseudogap state is a symmetric phase 
precursor state in the broken symmetry phase. In partic- 
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ular, the transition between the broken symmetry phase 
and the pseudogap state is not a phase transition but 
a cross-over prelude to the fully symmetric state, that 
the system enters when the lowest energy state of the 
effective potential is one where the modulus vanishes. 

It is reminded that the general arguments due to 
Kadanoff and Wilson imply, that in the ther¬ 

modynamical limit the phase transition properties of a 
material system are commonly universal, i.e. indepen¬ 
dent of the atomic level details. From this perspective, 
the construction of the phase diagram of a linear and chi¬ 
ral structureless homopolymer presented here, should be 
relevant for the understanding of the phase diagram of 
more elaborate linear chiral homopolymers, maybe even 
that of certain heteropolymers Indeed, linear poly¬ 

mers are presumed to have a very similar phase structure, 
quite independently of their chemical composition mi 
even though the phase where a particular polymer resides 
depends on many factors such as concentration, the qual¬ 
ity of solvent, ambient temperature and pressure. 

The article is organised as follows: The next Section 
describes the background and the methods. The stan¬ 
dard phase structure of linear polymers is first reviewed. 
The geometrical order parameter variables that are used 
to model the free energy of a homopolymer are defined, 
followed by an outline how the ensuing universal Hamil¬ 
tonian emerges in the limit of low wavelength deforma¬ 
tions. The order parameter that detects the presence of 
a pseudogap state is then defined. The zero temperature 
ground state is identified in the case of pure steric re¬ 
pulsion, and a universal attractive long-range interaction 
is introduced to model the effect of hydrophobic forces 
in the homopolymer chain. A detailed analysis of vari¬ 
ous Monte Carlo algorithms is presented, to identify one 
that is computationally most effective in the case of a 
single homopolymer chain. In the subsequent section 
the results are then described. The effect of the vari¬ 
ous parameters to the phase structure is revealed, and 
in particular the pseudogap state is identified. Finally, 
the phase diagram is constructed as a function of the 
various parameters. It is found that in the case of a chi¬ 
ral homopolymer the phase diagram has a much richer 
structure than in the case of a non-chiral homopolymer. 


II. METHODS 
A. Phases 

A review of the known phase structure of linear, non- 
chiral homopolymers is now presented, as a background 
and motivation for the subsequent study. 

Three different, universal phases are commonly iden¬ 
tified, and these phases are categorised by the way how 
the polymer structure fills the space mi- Under poor 
solvent conditions or at low temperatures, when the at¬ 
tractive interactions between the monomers dominate, a 
single polymer chain is presumed to collapse into a con¬ 


figuration which is space filling. On the other hand, in a 
good solvent or at high temperatures, when the repulsive 
interactions dominate and cause the chain to effectively 
swell, its geometric structure resembles that of a self¬ 
avoiding random walk (SAW). Between the two, there 
is a 0-regime (possibly a tri-critical 0-point) where the 
attractive and repulsive interactions cancel each other. 
In the 0-regime the polymer chain is presumed to have 
the characteristics of an ordinary random walk (RW). 
Finally, some polymers such as collagen for example, are 
more like straight, rigid rods. Each of these four phases 
- rigid rod, SAW, RW and the space filling one - can be 
characterised by the inverse of the Hausdorff dimension 
of the chain, called the scaling exponent [SHl]. This 
quantity is defined by the radius of gyration 


Rgyr — 




( 2 ) 


where the are the coordinates of the individual 
monomers. When the number of monomers N becomes 
very large, the radius of gyration has the asymptotic ex¬ 
pansion [24H261 [28] 


Rlyr RlN^^{l + RiN-^^ + ...) r-. (3) 

Here the length scale Rq, the Kuhn length, is the effective 
distance between the monomers in the large-limit. The 
Kuhn length is not a universal quantity, its value can in 
principle be computed from the atomic level details of the 
polymer and environment including pressure, tempera¬ 
ture and chemical microstructure of the solvent. The di¬ 
mensionless scaling exponent v i. e. the inverse Hausdorff 
dimension that governs the large-asymptotic form of 
equation (§, is presumed to be a universal quantity. Its 
numerical value is independent of the local atomic level 
structure of the polymer 0 BIMM]- The etc. are 
critical exponents and the Ri etc. are the corresponding 
amplitudes, and together they constitute the finite-size 
corrections. The etc. are universal quantities [22], but 
the i?i etc. are not universal [26] . 

The following mean field values are conventionally as¬ 
signed to V m- 


1/3 

collapsed 

1/2 

RW 

3/5 

SAW 

1 

rod 


Under poor solvent conditions or at low temperatures, 
the polymer collapses into a space-filling conformation 
[22IIHT] with the mean field exponent ly = 1/3. Folded 
proteins are commonly found in this phase. For an ordi¬ 
nary random-walk (RW) the mean field value is :/ = 1/2. 
This corresponds to the 6 regime, that separates the 
collapsed phase from the high-temperature self-avoiding 
random walk phase for which the Flory value ^' = 3/5 is 
found. Finally, when v = 1, the polymer loses its inher¬ 
ently fractal structure and behaves like a straight rod. 
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The transition between the collapsed phase and the 
SAW phase has been studied extensively, and it involves 
the RW phase as a tri-critical 0-point or more generally 
as a transitional 0-regime. But the transitions between 
the rigid rod phase and the other three, is less studied. 
However, there is a physically and biologically very im¬ 
portant scenario where such a transition could have a 
role; that of cold denaturation of a protein chain. The 
presence of all four phases Q opens the possibility of a 
4-critical point, under proper conditions [32]. 

Finally, the 0-point value = 1/2 is exact for a poly¬ 
mer with no long range interactions [^ . For a space-hlling 
structure the value = 1/3 is also exact, and similarly 

= 1 is exact for a straight, linear rod-like structure. But 
in the case of SAW the mean field value is corrected by 
fluctuations. A numerical Monte Carlo evaluation, com¬ 
puted directly by using the self-avoiding random-walk 
model on a square lattice, gives the estimate |26| 

V = 0.5877 ± 0.0006. (5) 

In the sequel the scaling exponent v in @ is evaluated 
as a function of various parameters in different phases 
of the homopolymer, using numerical simulations in the 
context of a universal off-lattice energy function. Of par¬ 
ticular interest is the effect of the parameter that charac¬ 
terises the chirality, and the parameter that characterises 
the strength of the attractive (hydrophobic) forces. 


B. Geometry 


The order parameters that determine the free energy 
of the homopolymer chain in terms of its geometry are 
now identified, following m- For this a homopolymer 
chain with i = monomers is considered, with 

the three dimensional space coordinates. The unit tan¬ 
gent vectors along the lines that connect two consecutive 
monomers are 


_ r,+i - r, 

|ri+i-ri| 

The unit binormal vectors are defined by 

1 ^2—1 ^ 

Di = ^-T 

|ti-l X til 

The unit normal vectors are defined by 


( 6 ) 

( 7 ) 


Hi = bj X ti (8) 

The three vectors (ni,bi,ti) determine an orthonormal 
frame at the monomer position ri. The discrete bond 
angles are 

Ki = K^+i^i = arccos (ti+i • ti) (9) 

and the discrete torsion angles are 


T, = Tj+i,i = sgn[(bj_i X bi) • ti] X arccos (bi+i • b*) 

( 10 ) 


Conversely, when the angles (fti, Ti) are known the dis¬ 
crete Frenet equation [27] 



/ cos hi cos r 

cos K sin r 

— sin hi\ 

i j 

(nA 

— sin r 

COST 

0 

• 

hA 

\ sin/s)COST 

sin hi sin r 

COSH ^ 


{uj 


( 11 ) 

determines the frames iteratively, by computing the 
frame at the position of the {i + 1)*^ monomer from 
the frame at the position of the monomer. Once all 
the frames have been constructed, the entire chain is ob¬ 
tained as follows. 


fc-i 


rfc = y] - Til • ti 


( 12 ) 


z=0 


With no loss of generality one can set Tq = 0, and orient 
to to point into the direction of the positive z-axis. 

A framing is necessary for the construction of the chain 
from the bond and torsion angles. But the equation (12) 


does not involve the vectors rii and bi. Thus any lin¬ 
ear combination of these two vectors could be chosen to 
define a framing, to construct the chain from the an¬ 
gles. This also determines the symmetry that enables 
the identification of the pertinent order parameter Q: 
Consider a local SO(2) transformation that rotates the 
frame (ni,bi) by an angle Ai leaving ti intact. 


( n\ / cos Ai sin Ai 0\ /n\ 

bj ( —sinAi cosAi oj- (bj (13) 

On the Frenet frame bond and torsion angles in ( [II] ), this 
has the following effect: 

(14) 


—>■ Ti ± Ai_i — A 


(15) 


where the {T°‘)bc = f-abc are the SO(3) generators 
[T“,T^] = eabcT‘^- The range of Ti is [—7r,7r) mod(27r). 
The equations (14) and (15) may be used to extend the 
range of the bond angle from [0, tt) to Ki into [—tt, tt) 
mod(27r). The extension is compensated for by the fol¬ 
lowing discrete Z 2 symmetry 


Kk ^ — Kk for all k > i 

Ti ^ Ti — t: 


(16) 


that leaves the chain intact. 

In the numerical simulations presented here, all the 
distances between nearest neighbour monomers are hxed 
to the uniform constant value 

|r,+i-r,| = <5 = 3.8 A (17) 

This equals the average distance between two consecu¬ 
tive Ca atoms along a protein backbone, measured in 
Angstrom’s. 
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A polymer is subject to steric constraints, due to over¬ 
lapping electron clouds and various short range Born 
repulsions. Accordingly the following forbidden volume 
constraint is introduced, 

|ri —= 3.8 A for \i — k\>2 (18) 


This is in line with the minimum distance observed be¬ 
tween any two Ca atoms, in folded protein structures. 
The numerical values (17) and (18) can both be indepen¬ 
dently modified, with no effect to conclusions. 

In the sequel only scaled dimensionless units are used, 
and in particular the dimensionless unit of length is one 
Angstrom. 


C. Free energy in the infrared limit 

The bond and torsion angles constitute a complete set 
of geometric variables, to describe protein Ca backbones 
[HH] . Furthermore, according to Eq. ^ and Eq. Q the 
structure dependent phase diagram of a homopolymer 
is determined by the three dimensional chain geometry. 
Thus the bond and torsion angles are a complete set of 
order parameters^ in the sense of Kadanoff and Wilson 
[ISHZS]: The geometrically defined, structural phase di¬ 
agram of a homopolymer can be fully determined by a 
thermodynamical free energy which is constructed from 
these order parameters only. 

A detailed derivation of the free energy used here is 
now presented. For this a homopolymer chain in thermal 
equilibrium is considered. Let F be the ensuing thermo¬ 
dynamical Helmholtz free energy. Thus, the minimum 
of F describes the chain configuration, under thermody¬ 
namical equilibrium conditions. The free energy is the 
sum of the internal energy U and the entropy 5, at tem¬ 
perature T 


F = U-TS (19) 

It is a function of all the inter-atomic distances 

F = Firafi) ; = |r„ - r,g| (20) 

where the indices a,/3,... extend over all the atoms in 
the homopolymer system, including those of the solvent 
environment. Consider the infrared, long distance limit 
where the characteristic length scales of spatial deforma¬ 
tions along the homopolymer chain around its thermal 
equilibrium configuration are large in comparison to the 
distance 0 between neighboring monomers. This is 
synonymous to an assumption that there are no abrupt 
wrenches and buckles along the chain, that there are only 
gradual long wavelength bends which is the limit of adi¬ 
abatic deformations. The completeness of the bond and 
torsion angles to describe a protein structure [33j implies 
that, in order to determine the thermodynamical phase 
state of the homopolymer chain, it is sufficient to consider 


the response of all the distances between all the atoms to 
the variations in the bond and torsion angles only, 

rap = rap{K,T). (21) 

Here, and in the sequel, (k, r) denotes collectively all the 
variables K,i and r^. 

Suppose that at a local extremum of the free energy, 
the bond and torsion angles along the homopolymer 
chain have the values 


(k*,a) = (Kio,Tio) (22) 


Consider a conformation where the {Ki,Ti) deviate from 
these extremum values. The deviations are 


AKi = Kj - KiO 
Ati — Tj T^o 


(23) 


Start by Taylor expanding the infrared limit Helmholtz 
free energy (19) around the extremum. 


F[rap=rap{ni,Ti)\ = F{k,t) = F{ko,to) 



dF 


dn 


Ank 


k |0 


dF 

dTk\ 


Ark 



1 d^F 

2 dnkdni 


Ak^Aki + 
\o 


d^F 

dKkTl |o 


AnkAriF 


1 n2 p 'I 

The first term in the expansion evaluates the free energy 
at the extremum. Since (KiOiAo) correspond to the ex¬ 
tremum, the second term vanishes. Denote in the sequel 
[Ki,Ti) collectively, as the variable pi. Then, 

F{k,t) = F{p) = 


1 (9^ F' 

= F{p^) + X ^ ApkApi + 0{A^) (25) 


k,l 


dpkdpi |o 


Following |34j the expansion (25) is re-arranged in terms 
of of the differences in the angles pi ^ r«), as follows: 

F{p) = ^{ Vk{pk-,Pok) + 


+ Zk{pk] P0k){pkPk+l + PkPk-l) + ...} (26) 

Here pofc denotes a combination of the various parameters 
iKio,Tio) along the chain. But Vkipk]Pok), Zk{pk]Pok) 
and so forth depend on the variable pk only on the site 
fc; these functions are ultralocal. The terms that are 
not shown explicitly, consist of higher order differences 
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PkPk+i with i > 2, and higher powers of the differences. 
The local terms Vk{pk) constitute the effective potential 

Veff=J2^k{Pk) (27) 

k 


The structure of the effective potential is commonly used 
to conclude whether a spontaneous symmetry breaking 
takes place [34]. 


The transition from (25) to (26) involves, a priori, an 


infinite re-arrangement of the terms in the Taylor ex¬ 
pansion (25). In particular, the expansion (26) has been 


designed so that in the continuum limit where distance 
between neighboring monomers vanishes i.e. d —> 0 in 
0. it becomes, at least naively, an expansion of the free 
energy in powers of momentum about the point where 
momentum vanishes: For a single scalar variable pk with 
continuum limit pk —>■ the corresponding continuum 

limit of (26) is the derivative expansion [31] 


FW = J 




(28) 


D. Effective Hamiltonian 


Clearly, the free energy must remain invariant under 
the local frame rotations (14), (15); the physical proper¬ 


ties of the chain do not depend on the choice of framing. 
Accordingly, it has been concluded [SHni that - in the 
unitary gauge - to the leading non-trivial order the free 
energy has the form 


N-l 

H = — 'y ) 2Ki^iKi -\- 
i=l 


i=N 

E 

2=1 


{2k? -k <7 (k? - nfff + cdK^rf} 


i=N 

E 

2=1 


-erf - QTi - bnjTi 


(29) 


This is adopted as the (effective) Hamiltonian, in the se¬ 
quel. In (29) q, m, a, b, c, d depend on the atomic level 


physical properties and the chemical microstructure of 
the homopolymer chain and its environment. In princi¬ 
ple, these parameters can be computed from this knowl¬ 
edge. (Note the combination cd in the last term of second 
sum; this choice is made for later convenience.) 

It can be shown |8HI3] that ( [2^ is the most general, 
universal and gauge i.e. frame rotation (13) invariant 


Hamiltonian, that models a homopolymer in the limit 
where the characteristic length scales of spatial deforma¬ 
tions around the minimum energy configuration become 
large in comparison to the distance 0 between con¬ 
secutive monomers. The effective Hamiltonian (29) co¬ 


incides with the naively discretized continuum Abelian 
Higgs Model Hamiltonian with one complex scalar field, 
when expressed in the unitary gauge and the U(I) gauge 
transformation is identified with the frame rotation (13); 


the term with parameter a is the Chern-Simons term, 
commonly introduced in gauge theories to break parity. 
In particular, (29) is unique in the sense of Kadanoff and 
Wilson. 


Implicit in (29) is the assumption that there are no 


abrupt wrenches and buckles along the polymer chain. 
Only small, gradual bends are present in the deviations 
around the energy minimum configuration, which is ob¬ 
tained by minimising the Hamiltonian (29). This defines 
the limit of adiabatic deformations. 

In line with the Abelian Higgs Model, see [3S] in the 


present context, the Hamiltonian (29) displays the dis¬ 


crete symmetry As in the case of the Abelian 

Higgs Model, this symmetry may become spontaneously 
broken by the ground state. It should be noted that (29) 
is not invariant under the local Z 2 gauge symmetry (16), 


as it coincides with the leading non-trivial contribution 
to an expansion of the Helmholtz free energy around a 
fixed background. To recuperate the Z 2 symmetry one 
may replace r in Hamiltonian by | sin 2r. Alternatively, 


the Hamiltonian (29) can be interpreted as a deformation 


of the standard energy function of the discrete nonlin¬ 
ear Schrodinger equation (DNLS) [MIIST]- The first two 
sums coincide with the energy of the standard DNLS 
equation, in terms of the discrete Hasimoto variable of 
[10] . The first (c) term in the third sum is the Proca 
mass that has a claim of gauge invariance; here the Proca 
mass is a “regulator”, as explained in [T3]. The second (a) 
term is the helicity, and the last (6) term is the conserved 
momentum. The last two terms break the Z 2 parity sym¬ 
metry, these two terms are responsible for helicity of the 
homopolymer chain. 

The simulations that are described in this article have 
been performed by keeping some of the parameter values 
fixed. In Table |l] the parameter values that are kept fixed, 
have been listed. The numerical values of these parame- 


TABLE I: The parameters in (291 that are kept fixed during 
our simulations. 


9 

m 

b 

c 

d 

3.5 

1.5 

0 

10“^ 

10"^ 


ters have been chosen in conformity with those, that are 
commonly encountered in the case of proteins [28l |38| ; 
for a protein, the torsion angles are much more flexible 
than bond angles. The value of m in Table [I] corresponds 
to an a-helical structure. 

The parameter a, which is not fixed, is of particular 
interest in the sequel. This is the parameter that breaks 
chirality; note that the momentum of the DNLS hierar¬ 
chy is not considered here i.e. b — 0 for simplicity. This 
term lacks a direct interpretation in the context of the 
Abelian Higgs Model. It turns out that the effects of 
this term are largely accounted for by the a dependent 
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helicity, in any case. 


Accordingly the minimum energy ground state of (29) is 


E. Pseudogap 

In [To] the following combination of the bond and tor¬ 
sion angles has been considered 


= tan ^ e*’’’ 


where the phase is 


N 




(30) 


(31) 


\k=l k=i+l 


The variable (301 is essentially the discrete version of the 


Hasimoto variable m, in terms of the Frenet frame coor¬ 
dinates: It is the complex variable that relates (29) into a 


generalised version of the discrete non-linear Schrddinger 
equation [36| |37|. It is also the present version of the 
complex order parameter 0 - Thus the pseudogap state 
can be identified with a state where the bond angles are 
non-vanishing and ordered 


(Ki) = Ko ^0 


(32) 


for some site independent kq, while torsion angles are 
essentially randomly fluctuating to the effect that 


(e"‘) « 0. 


(33) 


Accordingly in the sequel the pseudogap state is detected 
by monitoring both Ki and Ti simultaneously. It should 
be noted that since the effective potential (27) is insen¬ 


sitive to the phase in (30), the pseudogap state can be 


difficult to detect in terms of the minima of the effective 
potential alone. A dynamical computation that engages 
fluctuations is needed, to detect the presence of the pseu¬ 
dogap. 


It should be kept in mind, that a relation such as (32) 


is commonly deduced by inspection of the effective po¬ 
tential. In the full theory there are always corrections, 
due to fluctuations. In particular, in the full theory, at a 
finite temperature, (32) never vanishes identically. The 


modulus of the order parameter (30) is a positive definite 


quantity and thus, due to fluctuations, it always acquires 
a non-vanishing value in the full theory, as also shown in 
the simulations presented here. 


Ki = ± 


Ti = 


m2 - grf 


2c dm^ + l 


±m 


(36) 


In the sequel the parameter m has the fixed value, given 
in Table |T| throughout. As a consequence the ground 
state is controlled by the ratio a/c, and in the sequel the 
phase structure is investigated by varying the ratio a/c 
within the range a/c S [0,47r]. 


It should be noted that the configuration (36) does not 


necessarily describe the minimum energy homopolymer: 
For some parameter values there can be a conflict be¬ 
tween the values of {Ki,Ti) given by (36) and the forbid¬ 
den volume constraint (18). A configuration (36) which 


satisfies the forbidden volume constraint is a helix. But if 
the forbidden volume constraint is not obeyed, the lowest 
energy ground state configuration is the one that min¬ 
imises (29), subject to the constraint (18). 


The range of the bond angle become extended to neg¬ 
ative values, by the Z 2 symmetry (16). The two ground 
states Ki = ztm have the same energy. In addition of 
these two ground states, there can also be local minima 
of (29) that have the profile of a kink |38l [39] i.e. a 


domain wall that interpolates between the two ground 
states Ki = ±m. The energy of a kink is higher than the 
energy of the ground state helix ( jM] ). Two kinks can an¬ 
nihilate each other, thus any pair of kinks can be removed 
by continuous deformations of the chain. A single kink 
can be translated, so that it becomes removed through 
the ends of the chain. However, on a discrete lattice 
the translation invariance is commonly broken, by the 
Peierls-Nabarro barrier [TDI143j . Thus, in general it costs 
(thermal) energy to translate a kink along the chain. 

The present scenario is different from the one that ap¬ 
pears in the case of kinks in folded proteins (38] [39] . 


There, the parameter values in (291 are different for dif¬ 
ferent super-secondary structures i.e. helix-loop-helix 
motifs. A folded protein is described by a heteropoly¬ 


mer generalisation of (29), and the ground state is not a 
straight helix such as (36). A short analysis of a simple 


heteropolymer is presented in the sequel, in sub-Section 

[mDj 


F. Zero temperature 

The zero temperature ground state of the Hamiltonian 


(29) is a solution to the equations of motion. 


1 


(34) 


* c dKf + 1 
Ki+i = 2Ki — Ki-i + 2q{Ki — m‘^)Ki + cdr^Ki (35) 


G. Attractive interaction 


The constraint (181 is a purely repulsive interaction. 


and in the case of a homopolymer it models forbidden 
volume constraints which have a short spatial range. It 
could be generalized to include an attractive component, 
with a short range but one that exceeds the extent of 
forbidden volume constraints. Accordingly, the follow¬ 
ing rudimental extension of (18) to model both short 


range forbidden volume constraints and attractions is in- 
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troduced. 


U{r) = 


+00 


0 < r < (5 


t/o{tanh(r — Rq) — 1} 6 < r < +c» 


(37) 


Here <5 is the radius of the self-avoiding condition (18). 


For r < S the forbidden volume condition (18) persists. 


But for r > S there is a short range attractive interaction 
with strength determined by the parameter Uq; when the 


parameter Uq vanishes (37) reduces to (18). In the sequel 
this parameter will be varied, jointly with the ratio a/c 
that characterises helicity. 

The parameter Rq which determines the range of 
the attractive interaction, shall have the following value 
Rq = 5.0 A throughout; this choice is in line with 
all-atom molecular dynamics simulations where any long 
range interactions between atoms are commonly cut off, 
sharply, beyond distances around 10 Angstrom or so. 

The attractive interaction can be given a physical in¬ 
terpretation, in terms of “hydrophobic” forces: In the 
case of e.g. a protein chain under physiological condi¬ 
tions, there is an effective attractive interaction between 
those amino acids which are considered “hydrophobic” 
if]. Thus in the presence of the attractive interaction 
(37) the energy function (29) models a chain made of 


“hydrophobic” residues, with “hydrophobicity” that de¬ 
pends on the value of Uq. 

It is noted that qualitatively, the present results have 
been found to be quite insensitive to the details of the 


profile of the potential U{r). Accordingly (37) can be 
considered ’’universal”. 


H. Monte Carlo Algorithms 


The protein folding problem [46II48] is notorious for 
its computational complexity. A comprehensive all-atom 
simulation with Anton [ISl which is by far the fastest 
molecular dynamics machine available, can produce no 
more than a few micro seconds of a folding trajectory 
per day in silico, in the case of proteins with less than 
100 amino acids. Since many proteins take seconds, even 
days to fold into their native state starting from an initial 
random conformation, it could take thousands of years to 
fold such a protein with presently available computer re¬ 
sources. Even in the case of effective off-lattice models 


such as (29), (18) the simulation of a full folding trajec¬ 


tory is a formidable computational challenge, even with 
the most powerful computers available. Accordingly, to 
identify an effective computational scenario, the perfor¬ 
mance of three different Markovian Monte Carlo algo¬ 
rithms [44] have been tested. The aim has been to iden¬ 
tify a Monte Carlo algorithm that has the fastest rate 
of convergence towards a thermal equilibrium state, in 
the case of a single polymer chain. In these tests only 
the forbidden volume constraint (18) has been used, the 


effect of the attractive (“hydrophobic”) interaction has 
not been included. The three algorithms are 


1) Heat Bath algorithm 

2) Metropolis algorithm 

3) Mixed algorithm 

The curvature and torsion angles are updated accord¬ 
ing to a probability distribution, that satisfies the de¬ 
tailed balance condition 

P{{ ^new!t ^new} 7 { Hold , 'Told }) exp ( -/3i7 ( { Kold, T-old} ) ) 


— T’('[^T:old 1 "Fild}5 {^new5 "^new}) exp( /377('[/vn0^, T^ew})) 

Here /3 is the inverse Monte Carlo temperature. The 
equilibrium distribution 

exjp{—PH{K,T)} (38) 


of a canonical ensemble is obtained in the limit of an infi¬ 
nite number of updates. Each update consists of a “walk” 
through the entire chain with a provisional revision of 
each value {Ki,Ti) which is subject to the requirement 
that the forbidden volume constraint ( |18[ ) is preserved ; 
the three algorithms differ from each other only in the 
manner how the new values t”®'") are generated. 

It should be kept in mind in the sequel, that the Monte 
Carlo temperature T = j3~^ is not equal to the physi¬ 
cal temperature factor ksd where ks is the Boltzmann 
constant and the temperature 9 is measured in Kelvin 
scale. T is dimensionless quantity like energy H (see eq. 
(29) ) and all parameters in it (see Table |T]). Instead, in 
the low temperature collapsed regime general renormal¬ 
isation group arguments [IS] propose that dimensionless 
Monte Carlo temperature T is connected with real phys¬ 
ical temperature in the following way: 


InT = ksO -f ... 


(39) 


1. Heat Bath algorithm 

In the Heat Bath algorithm, new values (k"®"", t"®'^) 
are generated randomly, according to probability distri¬ 
butions 


and 


P{^Tn = — exp{-/3i7,,,«®-)} (40) 

^i.K. 


PirD = (41) 


Here and Hi r are the sum of all those terms in the 
Hamiltonian (29) that contain Ki and Tj, respectively, 
with the given index i. The and Zi^r are normali¬ 
sation factors. The updated values of k”®"' and rf®'^ do 
not depend on the previous values of Ki and r^. 

The probability density for k"®"" has the form 


P{Ki) ~ exp{-CiK| - C2Ki - C3 kJ 


(42) 














where 


In these tests, the following parameter values have been 
used in the hamiltonian (29), 


Cl = Pq 

C 2 = /3(2 - 2qm^ + - abn) (43) 

C3 = l3{-2{Ki+i + 


Thus (42) is non-Gaussian. On the other hand, the prob¬ 
ability density has the Gaussian profile 


Pin) ~ exp{-/?(^[d«;2 -h - ar,)} (44) 


Rejection sampling has been used to generate random 
numbers according to these probability distribution: Af¬ 
ter generating nf™ and rf®'" the forbidden volume con¬ 
dition (18) is checked, and the update is rejected when 
the condition is violated. 


2. Metropolis algorithm 

New values of k and r are generated according to Gaus¬ 
sian probability distributions, which is centered at the old 
values. The dispersion of each Gaussian can be adjusted, 
to enhance the convergence of the algorithm. The new 
values of k and t are accepted or rejected, in the same 
manner as in the conventional Metropolis algorithm. For 
example, in the case of Ti the probability of acceptance 
of a new value is 


P{Ti) = minjl, exp(—,0Ai/)} 


(45) 


where Ai7 is the difference of the energy between the 
new and the old configurations. In addition, the self¬ 
avoidance condition (18) is also verified at each step. 


3. Mixed algorithm 

The values of k"®'^ are generated in the same manner as 
in the Heat Bath algorithm, while for rf®'" the Metropolis 
algorithm is used. The convergence of the algorithm can 
be adjusted, by changing the dispersion of the Gaussian 
distribution in the r update. 


4- Algorithm comparison 

A priori, each of the three Monte Carlo algorithms 
should converge towards the same equilibrium distribu¬ 
tion albeit at a different speed. Thus, when the Marko¬ 
vian length is not sufficient and the equilibrium distribu¬ 
tion is not yet reached, the result depends on the algo¬ 
rithm and the number of Monte Carlo steps. Accordingly, 
the three algorithms have been tested and compared, to 
identify the one with the fastest convergence rate towards 
a known equilibrium state; for this, the zero temperature 
ground state described in sub-section II. E. is utilised. 


a = -1.0 X lO"'^ 
c = 1.0 X 10-4 


(46) 


These values have been chosen to reproduce a 
monotonous a-helical structure in a manner which is 


consistent with the forbidden volume constraint (18), as 


the lowest energy conformation; the choice is not unique. 
The simulations have been performed with varying chain 
lengths, from N=100 to N=900. In each case, simulated 
annealing has been used and the initial configuration is 
always a linear straight rod with Hi = Ti = 0. The ini¬ 
tial configuration is first heated to very high temperature 
values (up to T^ax = 100), where the structure is fully 
thermally randomised. This is followed by a slow cool¬ 
ing period, to the target temperature. The cooling takes 
place with small temperature steps, with each step equal 
to AT = 0.05 — 0.5 in logarithmic scale. After each step, 
10^ — 5 X 10^ Monte Garlo updates are performed along 
the whole chain, to ensure that it becomes thermalized 
to the ambient temperature. Here the term “termalized” 
means that the Marcovian chain reaches its equilibrium 
distribution. We will call “thermalization length” the 
number of Monte Garlo updates per one step of cooling 
process. The longer the chain, the longer the thermaliza¬ 
tion. For each temperature step only the last configura¬ 
tion is used for calculation of observables. Thus during 
one cooling process one final (thermalized) chain config¬ 
uration is obtained for each temperature value. The final 
phase diagrams have been calculated using 5000 Monte 
Carlo updates per one step of cooling, equal to AT = 0.05 
in logarithmic scale. One single cooling procedure takes 
6 CPU-hours, producing one configuration for each tem¬ 
perature value. Final statistics is compiled from 128 ther¬ 
malized chain configurations, for each set of parameters 
values. 


In the sequel the radius of gyration ([^ is utilised as 
the principal observable, to characterise the geometry of 
the homopolymer chain. Its value is calculated as an av¬ 
erage over statistically independent chain configurations, 
produced during a Monte-Carlo process. In Figure the 
value of the radius of gyration ([^ is compared, as a func¬ 
tion of Monte Carlo temperature for a chain with N=100 
monomers. The results obtained using the Metropolis 
algorithm are found to be very different from those ob¬ 
tained using the Heat Bath and Mixed algorithms; for 
small dispersion At the Mixed algorithm coincides with 
the Metropolis algorithm. But when At increases, the 
Mixed algorithm approaches the Heat Bath algorithm as 
shown in the Figure. 

There is no a priori reason why the results for the three 
algorithms should be different: The stationary distribu¬ 
tion is the same, in each of the three algorithms. But 
it is found that the Metropolis algorithm converges very 
slowly towards the equilibrium distribution. For this, the 
dependence of the result on the length of simulation has 
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FIG. 1: (Color online) Dependence of Rgyr on temperature 
for identical homopolymer chains simulated with different 
algorithms. Here At is dispersion of Gaussian distribution 
for generation of Tnew in mixed algorithm. In conventional 
Metropolis algorithm the dispersions are Ak — At = 0.01. 
Length of polymer chain is A^lOO, parameters of the Hamil¬ 
tonian are the ones from Table 1. 
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FIG. 3: Dependence of Rgyr on the thermalization length in 
Mixed algorithm. The number of updates, per one step in 
simulated annealing process (thermalization length), is shown 
along the horizontal axis. Parameters of this run are the same 
as for the figure 2. 




Thermalization length 

FIG. 2: Dependence of Rgyr on the thermalization length in 
Metropolis algorithm. The number of updates, per one step in 
simulated annealing process (thermalization length), is shown 
along the horizontal axis. Temperature is T = 10“^^, length 
of polymer A=100, and parameters of Hamiltonian are taken 
from the Table I. 


been analysed. In the case of the Metropolis algorithm, 
the results are shown in the Figure for T = 10“^^ 
which is the lowest Monte Carlo temperature value that 
has been used in the present simulations. As shown in the 
Figure, Rgyr continues to increase with increasing length 
of thermalization. The Metropolis algorithm approaches 
the Heat Bath algorithm very slowly, as a function of the 
simulation time and even at the present, relatively long 
simulation times the chain is still far from equilibrium 
distribution. 

On the other hand, the results shown in Figure 


FIG. 4: (Color online) Compactness index u at different tem¬ 
peratures for a homopolymer chain. Simulation performed 
with Heat Bath algorithm. Parameters of Hamiltonian are 
taken from the Table I. 


demonstrate that in case of the mixed algorithm the ra¬ 
dius of gyration Rgyr approaches a fixed value with in¬ 
creasing of the thermalization length. It is concluded 
that using the present simulation times, the Markovian 
homopolymer reaches a stationary distribution. Either 
the Heat Bath algorithm, or alternatively the Mixed al¬ 
gorithm with sufficiently large At, should be used to try 
and describe the thermal equilibrium configurations. 

The compactness index v i.e. the inverse of the Haus- 
dorff dimension has also been inspected, using the three 
different algorithms. The results for the Heat Bath al¬ 
gorithm are shown in Figure Between T = 100 and 
T = 1000 the value of v is essentially temperature inde¬ 
pendent, and apparently corresponds to the SAW phase. 
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At very low temperatures a transition to the rigid rod 
state, with z/ « 1, is observed. This result is in line with 
the general arguments that are presented in sub-section 
|IIG[ on the expected phase structure of the homopoly¬ 
mer model (29), (18). 


5. The final algorithm 


On the basis of the results from the test runs, the fol¬ 
lowing improved Heat Bath algorithm is employed in the 
sequel: The final probability distribution is 


1 


P = - exp{-/3[H + ^ {/(f, - fj)]} 


(47) 


i<j 


Here H is the Hamiltonian (29) and U is the potential 


(37). The Metropolis algorithm is used for acceptance. 


but with a proposal distribution that coincides with the 
Heat Bath algorithm: The new values of Ki and are 
generated using the distributions (|40|) and (ETl). The en¬ 


suing homopolymer configuration is then accepted, pro¬ 
vided it satisfies both the self-avoidance condition (18) 


and the Metropolis accept-reject condition that utilises 
the residual energy 


Eu = '^U(Ti-rj) 

i<j 


The acceptance criterion is 


exp{—/3AEu) > A 


(48) 


(49) 


where A is a random number which is uniformly dis¬ 
tributed between 0 and 1, and AEjj is the change in 
Ejj under the update of Ki and r^. 

Finally, the algorithm has been calibrated by consid¬ 
ering the limit of a truncated Hamiltonian , where the 
Hamiltonian El is removed and only the attractive poten¬ 
tial (37) together with the forbidden volume constraint 


(18) are retained. The numerical value of Uq determines 


solely the scale of the Monte Carlo temperature T, in 
the present simulations the value Uq = 15 is used. The 
results are shown in Figure A smooth transition is 
observed in the radius of gyration, from larger values at 
high temperatures to smaller values at low temperatures. 

The Figure|^shows how Rgyr depends on the length of 
polymer chain iV, and fitted to the leading order contri¬ 
bution in ([^. The results demonstrate how the equilib¬ 
rium distribution is reached with a sufficiently long ther- 
malization length, even in the region of low temperatures 
where the convergence is at its lowest. In particular it 
is found that the value of the compactness index n con¬ 
verges towards the mean field value 1/3 of the collapsed 
phase, as shown in Figure the difference between the 
space filling i/ = 1/3 and the numerically deduced v « 
0.36 is attributed to the hnite size corrections in ([^; the 
available computer power does not enable an identifica¬ 
tion of the amplitudes or the critical exponents 

i5i,... in ([^ . 
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FIG. 5: (Color online) Dependence of Rgyr on temperature T 
for a chain with 100 monomers, in the truncated model. Here 
n is the number of Monte Carlo updates per one simulated 
annealing step in temperature (thermalization length). The 
final algorithm described in paragraph H5 of the section II is 
used. 
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FIG. 6: (Color online) Dependence of Rgyr on the length of 
polymer chain N for different thermalization lengths with T = 
1, in the truncated model. The final algorithm is used. 


III. RESULTS 


The present variant of the Heat Bath algorithm has 
been used in extensive numerical simulations to investi¬ 
gate the phase structure of the homopolymer model 
(37). The results are summarised in Figures [8p^ 
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FIG. 7: Dependence of compactness index v on thermalization 
length when T = 1, in the truncated model. 

A. Effect of parameters 




logioT 


FIG. 8: (Color online) Dependence of Rgyr on temperature for 
a chain with 50 monomers. Parameters of the Hamiltonian are 
taken from the Table I. 


1. Parameter Uo 


Figures describe the properties of the radius of 
gyration for three representative values of the strength 
parameter Uq in (371, 


Un = 



(50) 


For the other parameters the values (46) are used. For 


each value of Uq, three different characteristic regimes are 
observed. In the high temperature limit the homopoly¬ 
mer is found in the SAW phase; this is confirmed in Fig¬ 
ure When the temperature decreases, the homopoly¬ 
mer enters a regime of decreasing Rgyr- Finally, there is 
the low temperature regime where the radius of gyration 


Rgyr has a small value; see Figure [T0| The compactness 


index v shows that when the thermalization increases the 
value of i' converges towards the values close to v ~ 0.39 
which is indicative of the mean field value 1/3. This 
is shown in Figure |11[ Again, the difference between 
the mean field value v = 1/3 and the measured value 
V « 0.36 is allocated to finite size corrections in @ ; the 
available computer power is not sufficient to deduce the 
detailed form of these corrections. 
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FIG. 9: Dependence of Rgyr on the length of chain for high 
temperature T = 10^'^ and with Uo = 0.01. Parameters of the 
Hamiltonian are taken from the Table I. 


v = 0.597 ±0.019 


the effect of such a change in a and c at finite temper¬ 
ature, in particular how it conspires with the parameter 
t/g, simulations have first been performed with 

a = c = 10-2 & Uo = 10-'‘ (51) 


It is concluded that when the scale for Uq exceeds that 
of a and c the low temperature state is a collapsed con¬ 
figuration. 


2. Parameters a and c 


According to (36) the classical ground state profile of 


remains intact when the parameters a and c are changed 
in such a manner that the ratio a/c is constant. To study 


with the values (46) for the remaining parameters. Thus, 
unlike in the previous simulations now the characteristic 
scale of the attractive interaction is smaller than that of 
the torsion angle dependent terms in the Hamiltonian. 

The results for the radius of gyration are presented in 
the Figure At high temperatures the chain is again in 
the SAW phase. Then, as temperature decreases, there 
is a transition to a regime akin the intermediate regime 
shown in Figure]^ Finally, there is a low temperature 
regime where the chain fluctuates around the classical 
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FIG. 10: (Color online) Dependence of Rgyr on the length of 
chain for low temperature T = 10“® and various thermaliza- 
tion lengths, and with Uq = 0.01. Parameters of the Hamilto¬ 
nian are taken from the Table I. 
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FIG. 11: Dependence of compactness index v on thermaliza¬ 
tion length. Temperature is equal to T = 10“® and Uq = 0.01. 
Parameters of the Hamiltonian are taken from the Table I. 


solution (36). The scale of transition to the low tem¬ 
perature regime is controlled by parameters a and c in 
Hamiltonian. 

It is concluded that when the scale Uq of the short- 
range attractive interaction is smaller than the scale of 
the parameters a and c, the low temperature limit is de¬ 
scribed by helical structures. 


B. Analysis of different phase regimes 

The bond and torsion angles form the complete set of 
local order parameters to probe the phase structure (|^, 
in the case of the present homopolymer model. These 
order parameters have the following characteristics, in 
the different regimes that have been analysed in Figures 


FIG. 12: Dependence of Rgyr on temperature. The chain has 


100 monomers, and parameters (a, c, Uq) are given in (511. 


•S 4 



•S 4 



FIG. 13: Distribution of (ni,Ti) in the high temperature 
SAW phase. Simulation parameters are the same as for the 
figure 8 {Uo = 0.01), temperature is T — 100. 


[8p^ the results are summarised in Figures [T3p6| 

In the very high temperature SAW phase, both the 
bond angle and the torsion angle are subject to large 
fluctuations; the simulation results are shown in Figure 
[T3| For the bond angles, the values are distributed in 


the range 0 < (k] < 


2 . 2 . 


The upper 
As 


limit reflects the forbidden volume constraint (18). 
temperature increases, the values of k become increas¬ 
ingly evenly distributed over this range so that in the 
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FIG. 14: Distribution of (ACi,ri) in the intermediate transi¬ 
tion regime: an example of pseudogap state (see its descrip¬ 
tion in section [ilE[ Simulation parameters are the same as 
for the figure 8 {Uo = 0.01), temperature is T = 1. 
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K T 


FIG. 15: Distribution of {Ki,Ti) in the collapsed phase. It 
shonld be observed that the distribution of Ti is asymmetric, 
corresponding to broken chirality. Simulation parameters 
are the same as for the figure 8 (Uo = 0.01), temperature is 
T = 10"'^. 



K T 


In the helical rod-like phase, both the bond and tor¬ 
sion angles become peaked around the classical values 
(36). The configurations are akin straight helical rods; 
a little like e.g. collagen when biologically active. The 
K distribution reflects the discrete Z 2 gauge symmetry. 
But the Z 2 symmetry in the r distribution observed in 
Figure is fully broken. 

The transition between the collapsed phase and the he¬ 
lical rod-like phase entails a transitition where the torsion 
angles become ordered. Due to the very low temperature 
values involved, the fluctuations are strongly suppressed 
and a simulation becomes tedious. It is conjectured that, 
when the temperature is kept in the low temperature 
regime, an initial helical rod-like structure but with pa¬ 
rameter values corresponding to the collapsed phase, is in 
a glassy phase. Vice versa, an initial collapsed configura¬ 
tion with parameter values in the helical rod-like phase, 
will eventually become subject to cold denaturation. 


C. Phase diagram 

The homopolymer phase is found to depend on three 
relevant scales. 


FIG. 16: Distribution of (Ki,ri) in the helical rod-like phase. 
Simulation parameters are the same as for the figure 12, 
temperature is T = 10“^. 


T —>■ 00 limit the distribution is fully uniform; the Fig¬ 


ure 13 shows the bond and torsion angle distribution at a 
generic but high temperature value. It is apparent from 
this Figure that both angles are disordered. 

It is concluded that the SAW phase is a disordered 
phase. 

Next, we observe the intermediate regime takes place 


-3 


< 


between temperature values within th e r ange 10 
T < 10^ as can be seen in FiguresandIn this inter¬ 
mediate region the values of the bond angle are found to 
become ordered. This is shown in Figure The values 
of Ki are thermally fluctuating at around |k| « 1.5. The 
values of the torsion angle remain largely disordered. The 
intermediate region is identified as a pseudogap state, as 
described in the Introduction. 

Finally, there are two different low temperature phases: 
The collapsed phase shown in Figure where the attrac¬ 
tive short-distance interaction dominates and the helical 
rod-like phase shown in Figure [T^ where the attractive 
short-distance interaction becomes weak. 

The distribution of bond and torsion angles in the col¬ 
lapsed phase are displayed in Figure The bond angle 
is highly ordered around the classical value (361 but the 
torsion angle remains disordered. However, there is an 
apparent spontaneous symmetry breaking that has taken 
place; the double well structure seen in the r distribution 
of Figure pH] has been removed, in a way that resembles 
the familiar spontaneous symmetry breaking in a Z 2 sym¬ 
metric potential well. 


- There is the extrinsic temperature scale where the 
values of Ki become ordered. This scale can also be con¬ 
trolled intrinsically, by the parameter q in (29), but the 
details have not been addressed here. 

- There is temperature scale where the values of Ti be¬ 
come ordered. This scale can be controlled intrinsically, 
by the parameter ratio a/cm (29). 

- Finally, the effects of the scale Uq for the short-range 
attractive interactions have been investigated. This 
parameter determines an intrinsic scale that controls 
the transition temperature alternatively to the collapsed 
phase, or to the helical rod-like phase. 


Thus, by changing the relations between the three 
scales the phase diagram of the homopolymer can be 
identified; the phase diagram is constructed here in terms 
of (T, a, Uq). All the remaining parameters are fixed, and 
given by the values in Table |l] 

The Figure shows the three-dimensional phase di¬ 
agram in the (T, a, C/q) space. The Figures ISpI show 
various cross-sections, taken at selected values of Uq i.e. 
these Figures show the phase diagram in Figure ini on 
the {a,T) plane, with different values of Uq. 


Figure [T^ shows the phase diagram, for a quite large 
value of Uq (strong coupling) At high temperatures there 
is the SAW phase. When temperature decreases there is 
the pseudogap state, that becomes the collapsed phase at 
low temperatures. Between the pseudogap state and the 
low temperature collapsed phase there is a d-regime, or 
rather a 0-point as it is observed only over a very narrow 
temperature range. 

- Figure [redisplays the phase diagram, when the value 
of Uq is lowered but still relatively large (intermediate but 
not weak coupling). At high temperatures there is again 
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FIG. 17: (Color online) The phase diagram on the 
(Uo, T, a) space. It was obtained for the polymer length 
N = 100 using our hnal algorithm. All parameters in 
the Hamiltonian except a are fixed according to the Ta¬ 
ble I. 


the SAW phase, followed by the pseudoogap state as the 
temperature decreases. At low temperatures the pseudo¬ 
gap state becomes converted either to the collapsed phase 
or to the straight rod phase, depending on the value of 
the helicity parameter a. In addition, there is a range of 
values of a, when a intermediate similar to the 0 -regime 
is observed between the pseudogap state and the straight 
rod phase. This is the 77 -regime. It is notable, that there 
is a possibility of a 4-critical point involving the pseudo¬ 
gap state, the 77 -regime, and the collapsed and straight 
rod phases. 

- Figurej^shows the phase diagram, as the value of Ug 
becomes further decreased (intermediate but not strong 
coupling). The collapsed state has entire disappeared, 
and replaced by the straight rod phase at very low tem¬ 
peratures. The 77 -regime displays a periodic structure 
in the parameter a. There appears to be a tri-critical 
point involving the pseudogap state, the 77 -regime and 
the straight rod phase. 

- Finally, in Figure 21 the weak coupling [/g phase 


diagram is displayed. The overall topology of the phase 
diagram is similar to the one in Figure [l 8 l but with the 
straight rod phase as the low temperature phase instead 
of the collapsed phase. 

The Figure!^ shows the phase diagram on the {Ug,T) 
plane, with helicity fixed. It is notable that there might 
be a 5-criticality involving the 77 - and 0-regimes, the pseu¬ 
dogap state and the collapsed and straight rod phases. 
However, the detailed investigation of this region of the 
phase diagram is beyond the capacity of the computer 
power which is presently available to us. 

In summary, it is remarkable how the phase diagram 
(figures FTpI ) is periodic in the helicity parameter a 


^ -2 
-4 


SARW 


0 ^•ogioUo=-3 

pseudogap 


o 


9-regime 


-6 

-8 


collapsed 


0.5 1.5 2.5 3.5 

a/(7i:c) 

FI G. 18: (Color online) A cross-section of the phase diagram 
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FIG. 22: (Color online) The phase diagram on the 
{Uo,T) plane, with a = 10“'*. 



logioT 


(when c is fixed). Moreover, there is a rapid transition 
between the helical phase where the compactness index 
ly = 1 and the collapsed phase where z/ « 1/3 at low tem¬ 
peratures, as shown in Figure There is a pseudogap 
state that appears as a transition regime, akin the con¬ 
ventional 61-regime, between the collapsed and the SAW 
phases. For the transition regime between SAW and heli¬ 
cal phases shown in Figure [20| there is a pseudogap state 
that should essentially coincide in its properties with the 
61-regime pseudogap state. But due to computational lim¬ 
itations the present analysis is not sufficient to confirm 
this. There is an apparent 4-critical, even 5-critical point 
as shown in Figures and but the detailed analysis 
of this region in the phase diagram needs to be performed 
using more extensive simulations, which is postponed to a 
future project. However, it is observed that the potential 
presence of a 4-critical point in a theoretical context very 
similar to the present one (z.e. Abelian Higgs model) has 
been reported in [5^ 


D. Heteropolymer and proteins 


Finally, it is inquired how the present results could 
be extended to heteropolymers, to draw conclusions on 
the potential phase structure of proteins. For this, the 
effect of perturbations that break the homogeneity of the 
homopolymer model have been investigated as follows: A 
collapse has been found to take place when the parameter 
Uq that characterises the strength of self-interaction is 
larger than the parameters a and c that characterise the 
torsion angle dependent terms. Thus, a short segment 
is introduced along the chain, where Uq is less than a 
and c. Accordingly, a simulation is performed where a 
heteropolymer is constructed so that for a short sub-chain 
of 12 monomers, the values of the parameters a and c is 
increased from a = c = 10“® to a = c = 10“^, while 
Uq = 10“* along the entire chain. Results of simulations 


are presented in Figures 23 and 24 for a chain with 150 
monomers. The dependence of Rgyr on temperature is 
shown in figure the phase diagram is very similar to 
For example, at low temperatures it is found that 


FIG. 23: (Color online) Dependence of Rgyr on temperature 
for a heteropolymer chain with 150 monomers Thermaliza- 
tion length n is number of updates per one step of simulated 
annealing. Parameters of the Hamiltonian and attraction po¬ 
tential are described in the paragraph D of the section III. 
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FIG. 24: Dependence of Rgyr on the length of polymer chain 
for low temperature T = 10“^ ®. Simulation parameters are 
the same as for the figure 23. 


v = 0.269 ±0.039 


the compactness index is close to the mean field value 
1 / = 1/3 of the collapsed phase. However, the geometry 
of a configuration in the collapsed phase is different: As 
shown in Figure a helical structure appears only in 
that sub-chain where the parameter values a and c have 
been increased. 


IV. SUMMARY 

The phase structure of chiral homopolymers under 
varying ambient temperature values have been investi¬ 
gated theoretically, in terms of a universal infrared limit 
energy function in combination with forbidden volume 
constraints (self-avoidance) and a short-range attractive 
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FIG. 25: (Color online) (left) Homopolymer does not dis¬ 
play helices in the ground state, (right) Addition of a sub¬ 
chain that breaks homogeneity, gives rise to a helix in the 
subchain region. Homopolymer configuration was obtained 
in simulation under the parameters inside the collapsed re¬ 
gion of the phase diagram (see hgure 17). In case of het¬ 
eropolymer the simulation parameters are the same as for 
the figure 23, dimensionless temperature T (see paragraph 
HH) is equal to 


interaction between residues. As such, the model should 
provide a realistic description even in the case of het¬ 
eropolymers that display an approximatively repeating 
monomer pattern, provided the scale of the repeat can 
be considered small in comparison to the chain length. 
A biologically important example is given by collagen, 
the most prevalent protein in a human body, in which 
case the short range attractive interaction models weak 
hydrophobicity of the amino acids. 

It is found that the phase diagram displays a high level 
of complexity, in terms of the parameters that control the 


helicity, and the strength of the attractive interaction. 
In particular, the low energy phase is either like a linear 
one dimensional straight rod, or a space filling collapsed 
configuration. At intermediate temperatures, there is a 
state which can be identified as an example of the pseudo¬ 
gap state, and there are also intermediates that are more 
like the conventional 0-regime. It is possible that these 
regimes merge, in the thermodynamical limit, at least for 
some range of parameter values. However, the possibility 
of the existence of 4-critical, even 5-critical points in the 
phase diagram is also proposed, but can not be confirmed 
with presently available computer power. 

The extension of the present approach to investigate 
properties of proteins and other heteropolymers remains 
a challenge to future research. 
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